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In the first paper of this series it was found that the jy-spin 1/2 holons, spin 1/2 spinons, and c pseu- 
doparticles whose occupancy configurations describe the energy eigenstates of the one-dimensional 
Hubbard model emerge from the electron - rotated-electron unitary transformation. An important 
breakthrough is that the theory provides relevant information about the relationship of the original 
electrons to these quantum objects. In this second paper we discuss and clarify how such a relation 
can be used in a program for evaluation of finite-energy few-electron spectral functions. As a first 
step, here we characterize the dominant holon and spinon microscopic physical processes that orig- 
inate more than 99% of the few-electron spectral weight. These dominant processes are related to 
exact selection rules for the values of the number of holons and spinons generated or annihilated by 
application onto a ground state of rotated-electron operators. While our theory also describes the 
higher-order microscopic processes associated with the remaining less than 1% few-electron spectral 
weight, the clarification and study of the above dominant processes is valuable and useful for the 
further understanding and description of the few-electron spectral properties observed in real low- 
dimensional materials. Moreover, in this paper generalize the concepts of a lower Hubbard band 
and upper Hubbard bands to all values of on-site Coulombian repulsion. 
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> ■ I. INTRODUCTION 

o\ 

' In this paper and its two companion papers HQ the relation of the electrons of the one-dimensional (ID) Hubbard 
, model to the quantum objects whose occupancy configurations describe its energy eigenstates is investigated for the 
whole Hilbert space. The study of such a non-perturbative relation is a necessary step for the description of the 
CO \ finite-energy few-electron spectral properties of the many-electron quantum problem. 

The study of the one- dimensional (ID) Hubbard model H S H S S S Q for interacting electrons is of general 
importance because the understanding of correlated systems and of their finite-energy excitations is still far from 
complete. The problem of the ID Hubbard model in the limit of large and infinite on-site Coulomb repulsion was 
previously studied in the literature by many authors |Il[l3[ll[l!H[llH[llHH Our studies of 

the model for finite values of the on-site repulsion U are motivated by the anomalous one-electron and two-electron 
spectral properties observed in metallic and insulating phases of quasi-one-dimensional materials, which cannot be 
described by the usual Fermi-liquid theor y |23u24l. Recently there has been a renewed exp erimental interest on the 
properties of these materials [Hi! El El IMiflll El El IM El El El El E3 IMIpiSlEiypJii . Some of these 
experimental studies observed unusual finite-energy/frequency spectral properties 28, 29, 30l l32ll41| . Since in the case 
of finite-excitation energy the Luttinger liquid description does not apply |4fjl l47l l48j] . these finite-energy/frequency 
spectral properties are far from being well understood. However, there are indications that electronic correlation 
effects might play an important role in the finite-energy physics of these low-dimensional materials [2I El El El El • 
For low values of the excitation energy the microscopic electronic properties of these materials are usually described 
by systems of coupled chains. On the other hand, for finite values of the excitation energy larger than the transfer 
integrals for electron hopping between the chains, ID lattice models taking into account the screened electron-electron 
Coulomb repulsion are expected to provide a good description of the physics of these materials. The simplest of these 
models is the above ID Hubbard model USSSji which describes such electron-electron interactions by an effective 
on-site Coulomb repulsion U. This model corresponds to a non-perturbative electronic problem. According to the 
results of the companion paper its energy eigenstates can be described by occupancy configurations of holons, 
spinons, and c pseudoparticles. Holons and spinons has also been studied for other models |43, |5(J, |5ll |52|- The 
ID Hubbard model is often considered a suitable model for the description of the electronic correlation effects and 
the non-perturbative microscopic mechanisms behind the unusual few-electron spectral properties observed in quasi- 
one-dimensional materials [28l I29I l4ll |44| . Moreover, recent angle-resolved ultraviolet photoemission spectroscopy 
revealed very similar spectral fingerprints from both high-T c superconductors and quasi-one-dimcnsional compounds 
[26| . The similarity of the ultraviolet data for these two different systems could be evidence of the occurrence of a 
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charge-spin separation associated with holons and spinons. The anomalous temperature dependence of the spectral 
function could also indicate a dimensional crossover |26t lo-'l |5J, |55| • The results of Refs. |5y, |57| also suggest that 
the unconventional spectral properties observed in two-dimensional (2D) materials could have a ID origin. Thus the 
holons and spinons could play an important role in spectral properties of both ID and 2D low-dimensional materials. 

In this second paper we continue the studies of the first paper of this series, Ref. 0] . As a preliminary application of 
the connection of the concept of rotated electron to the quantum numbers that label the energy eigenstates provided by 
the Bethe-ansatz solution and ry-spin and spin symmetries, in this paper we use exact holon and spinon selection rules 
for rotated-electron operators in the study of the holon and spinon contents of few-electron excitations. This reveals 
the dominant holon and spinon microscopic physical processes that generate more than 99% of the spectral weight 
of few-electron excitations. While our theory also describes the higher-order processes associated with the remaining 
less than 1% electronic spectral weight, the clarification of the dominant holon and spinon microscopic mechanisms 
is valuable for the further understanding and description of the few-electron spectral properties observed in real 
low-dimensional materials. Fortunately, a preliminary application of the theoretical predictions which follow from our 
study of the holon and spinon contents of few-electron excitations leads to quantitative agreement with the charge and 
spin spectral branch lines observed by means of angle-resolved photoelectron spectroscopy (ARPES) in the organic 
conductor TTF-TCNQ 58j]. The preliminary results reported in Ref. j5^ confirm that from the experimental point 
of view only the spectral weight associated with the dominant holon and spinon microscopic processes is observed. 
Moreover, in this paper we introduce the concept of an effective electronic lattice. The expression of the electrons in 
terms of holons, spinons, and pseudoparticles through the electron - rotated-electron transformation studied here for 
all values of U is a first necessary step for the evaluation of finite-energy few-electron spectral function expressions, 
as further discussed in Sec. V. 

In this paper we also discuss and clarify how the relationship of the original electrons to the quantum objects whose 
occupancy configurations describe the energy eigenstates can be used in a program for evaluation of finite-energy 
few-electron spectral functions. The successful fulfilment of such a program needs the concepts of local pseudoparticle 
and effective pseudoparticle lattice introduced in the third paper of this series, Ref. 0- Moreover, in this paper 
we generalize the concepts of a lower Hubbard band and upper Hubbard bands to all values of on-site Coulombian 
repulsion. 

The paper is organized as follows: In Sec. II we provide a short introduction to the ID Hubbard model. The 
definition of the upper Hubbard bands in terms of rotated-electron double occupation, the holon and spinon selection 
rules for rotated-electron operators, and the concept of an effective electronic lattice are presented and introduced 
in Sec. III. The dominant holon and spinon microscopic physical processes that control the few-electron spectral 
properties are studied in Sec. IV. Finally, in Sec. V we present the discussion of our results and the concluding 
remarks. This includes the discussion of the application of the concepts introduced in the present series of three 
papers to the fulfilment of a program for evaluation of finite-energy few-electron spectral functions. 

II. THE ID HUBBARD MODEL 

In a chemical potential /i and magnetic field H the ID Hubbard Hamiltonian can be written as, 

H = H H -(U/2)[N + N a /2}+ J2 (1) 

OL — C. s 

where, 

H H =f + UD, (2) 

is the "simple" Hubbard model. The operators 

N a 

f = -*££ £ ( 3) 

i=i <r=Ta5=-i,+i 

and 



N a 



(4) 
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on the right-hand side of Eq. J3J) are the kinetic energy and the electron double occupation operators, respectively. 
The operator 



ftj, a — c j,a c j, a , (5) 

of Eq. (J2J counts the number of spin a electrons at real-space lattice site j = 1, ...,N a . The number of lattice sites 
-ZV a is even and large and N a /2 is odd. We consider periodic boundary conditions. The operators cj a and Cj, „ which 
appear in the above equations are the spin a electron creation and annihilation operators at site j. Moreover, on the 
right-hand side of Eq. Q), /i c = 2/i, fi s = 2/jqH, hq is the Bohr magneton, and the operators S'i and S s z are given in 
Eq. (2) of Ref. and are the diagonal generators of the 77-spin S c and spin S s SU(2) algebras 0,IEIE3; respectively 
The electron number operator reads N — where the operator, 

3 

counts the number of spin a electrons. The momentum operator reads 



P 



N a 



-j+l, °Ji CT 



(7) 



and commutes with the Hamiltonian introduced in Eq. (|T|> . 

There are spin-up electrons and N± spin-down electrons in the chain of N a sites, lattice constant a, and length 
L = [N a a] associated with the model Q ■ We introduce Fermi momenta which in the present thermodynamic limit 
L — > 00 are given by ±/cf<t = ±7r?V and ±kp = ±[&ft + fcpjJ/2 = ±im/2, where n a — N a /L and n — N/L. The 
electronic density can be written as n — n-f + and the spin density is given by m — — ny. In general we consider 
densities in the domains < n < 1/a ; l/a < n < 2/ a and — n < m < n ; — (2/ a — n) < m < (2/a — n), respectively. 
However, in the case of the study of transitions whose initial state is a ground state, for simplicity we restrict our 
considerations to values of the electronic density n and spin density m such that < n < l/a and < m < n, 
respectively. The Hamiltonian Hso(4) = Hh — (U/2) [N + N a /2) commutes with the six generators of the 77 spin 
S c and spin S s algebras their off-diagonal generators being given in Eqs. (9) and (10) of Ref. 0. The 

Bethe-ansatz solvability of the ID Hubbard model is restricted to the Hilbert subspace spanned by the lowest-weight 
states (LWSs) of the 77-spin and spin algebras, i.e. such that S a = —S" 59]. 

As in the first paper of this series [ij , here we denote the holons (and spinons) according to their 77-spin projections 
±1/2 (spin projections ±1/2). For the definition of the holon and spinon numbers, their relations, and other aspects 
of the holon, spinon, and pseudoparticle description which are useful for the studies of this paper see the companion 
paper [jj- An important concept introduced in Ref. is that of CPHS ensemble subspace where CPHS stands for 
c pseudoparticle, holon, and spinon. This is a Hilbert subspace spanned by all states with fixed values for the —1/2 
Yang holon number £ Ci -i/2, —1/2 HL spinon number L S) _i/2, c pseudoparticle number N c , and for the sets of a, v 
pseudoparticle numbers {N c ,v} and {iV S)J/ } corresponding to the v = 1,2,3, ... branches. We note that according to 
the notation of Ref. in HL spinon HL stands for Heilmann and Lieb. 



III. THE UPPER HUBBARD BANDS, SELECTION RULES FOR ROTATED-ELECTRON 
OPERATORS, AND THE EFFECTIVE ELECTRONIC LATTICE 



In this section we generalize the large-t/ concept of upper Hubbard band [l2| to all values of the on-site electronic 
repulsion U in terms of rotated-electron double occupation. Moreover, we use the relation of rotated electrons to 
holons and spinons to introduce rotated-electron exact selection rules. These rules limit the values of the deviations 
in holon and spinon numbers associated with excited states generated by application of rotated-electron operators 
onto any eigenstate of the electron spin a number operator |JBJ. Deviation values outside the ranges given by these 
selection rules correspond to forbidden excited states. In addition, in this section we introduce the concept of an 
effective electronic lattice. This concept is further used in the third paper of this series, Ref. 0, in the introduction 
of the concepts of a local pseudoparticle and an effective pseudoparticle lattice. The latter concepts play an important 
role in the studies of finite-energy few-electron spectral functions of Refs. [60l l6l| . 
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A. ROTATED ELECTRONS AND THE UPPER HUBBARD BANDS 

As discussed in the companion paper Q|, the operator cj a represents the rotated electrons. The relation of such 
an operator to the electronic operator c] a is such that, 

cl a = V{U/t)cl a V\U/t). (8) 

Here V(U/t) is the electron - rotated-electron unitary operator uniquely defined by Eqs. (43)-(45) of Ref. 0. The 
rotated-electron double-occupation operator plays an important role in our studies and is given by, 



D = Vl (U/t) D V(U/t) = £ c 3 t T cj, T 4, x Cj, x , (9) 

3 

where D is the electron double occupation operator given in Eq. Q . Note that cj a and ct ct are only identical in the 
U/t — > oo limit where electron double occupation becomes a good quantum number. 

Let the operator Ojy(k) (or Ojy(k)) have momentum k and be a product of a finite number 



of one-electron (or one-rotated-electron) elementary creation and/or annihilation operators. Here the ratio j\f/N a 
vanishes in the present thermodynamic limit and Ni c j 3 is the number of electron (or rotated-electron) creation and 
annihilation operators for l c — — 1 and l c — +1, respectively, and with spin down and spin up for l s — — 1 and l s = +1, 
respectively. 

The studies we present below in this paper refer to the transformations generated by the application of the general 
TV-electron operator Ojs/{k) (or A/"-rotated-electron operator Oj\[(k)) onto any eigenstate of the spin a electron number 
operator ©. However, we focus our attention on the particular case of states Oj\f{k)\GS) (or Oj\r(k)\GS)) which 
result from application of such an operator onto a ground state \GS). We note that the number of spin a electrons 
equals the number of spin a rotated electrons. 

According to the results of the companion paper 0, rotated-electron double occupation equals the number of —1/2 
holons. We emphasize that for electronic densities n such that < n < 1/a the ground state has no —1/2 holons 
Since rotated-electron double occupation is a good quantum number, in this case we can write a general jV-electron 
spectral function as follows, 



oo 

4v(fc, w) = E Ai*M k > w ) ' ( n ) 

M=0 

where 



A AfM (k 1 = £ \{M,j\djsf{k)\GS)\ 2 S(u> - uud) . (12) 

3 

the excitation energy u> is positive, the M summation of Eq. 1)11(1 runs over the values of rotated-electron double 
occupation M — 0,1,2,... of the excited energy eigenstates \M,j), the j summation of Eq. fT^l runs over all 
available excited energy eigenstates \M,j) with the same value M of rotated-electron double occupation, and cJmj = 
[Em.] — Eqs] are the excitation energies relative to the initial ground state. Examples of few-electron operators Oj\f(k) 
whose spectral-function excitation energy is usually considered positive are the one-electron spin a creation operator 
of momentum k, the two-electron charge operator of momentum k, and the two-electron singlet or triplet Cooper pair 
addition operator of momentum k. 

We call lower Hubbard band the spectral weight distribution associated with the function Aj\ffi(k, uj) defined in 
Eq. (|12|) . Moreover, we call Mth upper Hubbard band the spectral weight distribution associated with the function 
Aj\f,Ai(k, u>) defined in Eq. (|12fl such that M > 0. In the ensuing section we find that for one-electron and two- 
electron operators most spectral weight correspond to the lower Hubbard band and first (or first and second) upper 
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Hubbard band(s). The detailed study of the few-electron spectral functions Ijlll) is for finite values of U a complex 
many-electron open problem. Such a problem is studied in Ref. [6l| by the combination of the pseudofermion 
representation introduced in Ref. [§D] with the holon and spinon description introduced and studied in this paper and 
in its companion papers 0,0] , as discussed in Sec. V. In this paper we use the relation of electrons and rotated electrons 
to the quantum objects whose occupancy configurations describe the ground state and excited states \M, j) that appear 
in Eqs. and i|12|) to obtain useful information about the finite-energy few-electron spectral distributions. 

An interesting property is that in the metallic phase of the ID Hubbard model the minimum value of the excitation 
energy uj for the lower Hubbard band is zero. This justifies the designation lower band. On the other hand, the 
minimum value of the excitation energy for the Mth upper Hubbard band is for finite values of U finite. Such a value 
depends on the minimum excitation-energy value for the creation of a —1/2 holon. In the case of the Hamiltonian 
(|T|). this minimum energy value is for electronic densities n such that < n < 1/a given by AE = E u . For a zero 
spin density m — initial ground state, the energy E u = 2/z (U /t, n) equals twice the chemical potential and is an 
increasing function of the on-site repulsion U with the following limiting values, 

E u = it cos(7rna/2) ; U/t^>-0, 

= U + 4t cos(Tma) ; U/t^oo. (13) 

Moreover, this energy parameter for any value of U/t is a decreasing function of the electronic density n such that, 

E u = U + At; n^O, 

= E MH ; n 1 , (14) 

where Emh is the half-filling Mott-Hubbard gap Q. On the other hand, if one uses as initial state the same m = 
ground state but considers the excitation energy of the "simple" Hubbard Hamiltonian (J2J), there is a shift in the 
definition of the ground-state zero-energy level. For this Hamiltonian the minimum excitation-energy value for creation 
of a —1/2 holon is simply given by U. 

Since rotated-electron double occupation equals the number of —1/2 holons, the minimum excitation-energy value 
for an excited energy eigenstate with rotated electron double occupation M' > is AE — M' E u . This justifies why 
all gapless excitation branches are associated with zero rotated-electron double occupation states. The excitation 
energy AE — M' E u is measured from the ground-state level and corresponds to the Hamiltonian . If instead 
we consider the ground-state zero-energy level associated with the Hamiltonian the value of such an excitation 
energy becomes M' U. 

Any operator O can be expanded as follows, 

0=J2 °ME* , (15) 
Me z 

where M = ... — 3, —2, —1, 0, 1, 2, 3, ... is such that application of the operator Ome u onto any eigenstate of rotated 
electron double occupation of eigenvalue M' leads to a final state with rotated electron double occupancy [M'+M] > 0. 
Also the corresponding rotated operator O such that, 

6 = v ] {u/t)6v(u/t), (16) 

can be decomposed as, 

6 = Ome u ■ (17) 

The index ME U of expressions (|15H and l|17|l plays the same role as the index nU of Ref. While M — n, the 

energy E u only equals U in the U/t — > oo limit. Since our study refers to all values of U ft, here we use E u instead of U 
to label the energy of the different upper Hubbard bands defined above. For finite values of U jt the choice of Ref. 0] 
refers to the zero-energy level of the Hamiltonian (0) , whereas our choice corresponds to the ground-state zero-energy 
level of the more suitable Hamiltonian Our choice is justified by the fact that in real experiments the excitation 
energy is, in general, measured from the Fermi level. Our notation tells that the minimal amount of excitation energy 



6 



for a transition which results from application of an operator component Ome u onto any eigenstate of rotated electron 
double occupation is ME U . This excitation energy can be positive (M > 0) or negative (M < 0). We note that 
application onto a state of rotated-electron double occupation M' < | M\ of a M < operator component Om_e„ 
gives zero. Since for the present case of electronic densities n such that < n < I /a, the ground state has zero 
rotated-electron double occupation, all the final excited states of the general ./V-electron spectral function l(ll|l - l(12|) 
have M > 0. 

The operators V'{U/t) and V(U/t) associated with the electron - rotated-electron unitary transformation can be 
written as, 



VHU/t) = ; V(U/t) = e s . (18) 

These operators are uniquely defined by Eqs. (44) and (45) of the companion paper p| and can be expressed in the 
form of an expansion in t/U for the operator S An important role in such an expansion is played by the rotated 
kinetic energy operator T, associated with the operator J2J). The decomposition ltT7|l of this operator leads to three 
terms, 



f = f + f Eu + f_ Eu , (19) 



where 



N. 

-t 

j=l cr=T,l<5=-l, + l 



Y Y Y [^-ni'^^eC-j+S^il-hj+S.a) 

(20) 



+ "'J, 5 Cj, a 5j+<5, cr flj+S, S 



te u = -t ^2 ^2 "•j. s 5 J, ° d J+ s ' ° ~ ™j+ s > s ) • 

J=l <r=T4 (5=-l,+l 



(21) 



and 



f_£; tt = -f^^ ^ (l-n 3 -, ff )ct CT c,' +< 5, CT nj +a , s . (22) 
j=i o-=t,1 *=-l,+l 

In these expressions a is such that | =J, and J, ="f. The above-mentioned expansion of the operator S of the expressions 
of Eq. I|18|) in powers of t/U is obtained by systematic use of Eqs. (44) and (45) of the companion paper pj, with 
the result, 



(23) 



where the components of the rotated kinetic energy operator are given in Eqs. (|20|l - 122|) . While equations (44) and 
(45) of Ref. uniquely define the electron - rotated-electron unitary operator V{U/t) of Eq. 118(1 for all values of 
U/t, it is difficult to extract a closed form expression for that operator from these equations. 

Let us consider the case of few-electron operators Otfik) whose spectral- function excitation energy is usually 
considered negative. Example of such operators are the one-electron spin a annihilation operator of momentum k and 
the two-electron singlet or triplet Cooper pair annihilation operator of momentum k. We note that according to the 
results of the companion paper 0, the number of rotated-electron empty sites equals the number M c +1 / 2 of +1/2 
holons. According to Eq. (C24) of the same paper, for values of the electronic density n such that < n < 1/a the 
ground state has a number of +1/2 holons whose value is given by M? +] y 2 = [N a — N]. Since the number of rotated- 
electron empty sites is a good quantum number, we can write a general A/"-electron spectral function associated with 
such a type of few-electron operators as follows, 



M=0,±l,±2,... 



(24) 
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where 

w) =E KM,j|0^(fc)|G5)| 2 «(o; + w A)J .) , (25) 

the excitation energy ui is negative, the M = [M c ,+1/2 — +1 / 2 ] summation of Eq. 1|24|) runs over the values of 
the deviation in the number of rotated-electron empty sites M = 0, ±1, ±2, ±00 of the excited energy eigenstates 
\M,j), the j summation of Eq. (|25|l runs over all available excited energy eigenstates \M,j) with the same value 
M Cj +1/2 = [M® + i/2 + M] f° r the number of rotated-electron empty sites, and Wjg-j = [£-jT/.j — -^gs] are the excitation 
energies relative to the initial ground state. Note that in the case of half filling the ground-state value for the number 
of rotated-electron empty sites is zero and the M summation of Eq. I|24|l runs over positive values of rotated-electron 
empty sites M — 0,1, 2, ...,00 only. Moreover, for electronic densities n > l/a the role of the numbers of rotated- 
electron doubly occupied and empty sites is interchanged. For these densities the summations over the values of 
rotated-electron doubly occupied and empty sites of Eqs. (ftl"]) and J23J, respectively, correspond to both positive and 
negative integers and to positive integers, respectively. This is because for such densities the ground state has no 
+1/2 holons. In this paper we consider mostly the case of electronic densities n < l/a, where the ground state has 
no —1/2 holons. 

As further discussed in Sec. V, the first step for the evaluation of the general few-electron spectral functions (f 1 1|> — 
HI 21 and (|24|I - H25|I is the expression of the TV-electron operator Oj^(k) in terms of the corresponding rotated operator 
CV(/c) defined in Eq. (JTBJ) as follows [6l|, 

6 u {k) = V{U/t)6 M {k)V\U/t) 

= 6^{k) + [S, 6 H {k) ] + \[S, [S, 6 M {k) ]] + .... (26) 

According to the results of the companion papers [J Q , the relation of the rotated operators on the right-hand side of 
Eq. (|26|l to the quantum objects whose occupancy configurations describe the ground state and excited states \M, j) 
and \M, j) of the few-electron spectral functions l|ll|) - (|12|) an( l l|24ll - 125|) . respectively, is well defined. Expression 
of the second expression of Eq. I|2t)|l in terms of the creation and annihilation operators for these quantum objects 
plays a crucial role in the evaluation of the matrix elements of the few-electron spectral functions l|12|) and l|25|l . This 
justifies the importance of the holon, spinon, and c pseudoparticle description introduced in the companion paper Q 
and further studied here and in the companion paper [2| for the study of the few-electron spectral properties at finite 
values of the excitation energy. Solution of this problem is a direct application of the relation of rotated electrons to 
the quantum numbers of the Bethe-ansatz solution and 77-spin and spin symmetries found in the companion paper 
0. Such a breakthrough provided the relation of electrons and rotated electrons to the above-mentioned quantum 
objects. 

In the ensuing section we find that the first operator term of the second expression on the right-hand side of 
Eq. (|26fl corresponds to the dominant microscopic physical processes that control the few-electron spectral-weight 
properties. In the case of JV = 1-electron operators we find that more than 99% of the electronic spectral weight 
generated by application of the operator Oj^ onto a ground state corresponds to application of the corresponding 
A/"-rotated-electron operator Ojy onto the same state. This latter operator has a simple expression in terms of holon 
and spinon elementary operators, and it is associated with the dominant holon and spinon microscopic mechanisms 
that control the A/"-electron spectral properties. 

B. EXACT HOLON, SPINON, AND c PSEUDOPARTICLE SELECTION RULES FOR 

ROTATED-ELECTRON OPERATORS 

Let \(j>) be any eigenstate of the spin a electron number operator © and Oj^ (or Oj^f) a general A/"-electron (or 
A/"-rotated-electron) operator. The corresponding state 0_\f\(p} (or Oj^\(f))) is also an eigenstate of the spin a electron 
number operator JBJ- Such a state belongs to an electron ensemble space whose electron numbers differ from the 
numbers of the initial state by deviations AiVj and AiV^ , such that AN = AN^ + AN± . Since the number of spin 
a electrons equals the number of spin a rotated electrons, these deviations also refer to the latter numbers. These 
deviations can be expressed in terms of the numbers A/i , i s of Eq. (|10fl as follows, 



A7V T = ]T (-UA/L+1 
; c =±i 



AN i = E 

; c =±i 



(27) 
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Following the results of Refs. PJH|, let us consider the four expectation values R a .i a = {R a .i a ), where a = c.s 
and l a = — 1, +1 and the corresponding operator R c ,—i counts the number of electron doubly-occupied sites, R c ,+i 
counts the number of electron empty sites, i? Si _i counts the number of spin-down electron singly-occupied sites, and 
R s ,+i counts the number of spin-up electron singly-occupied sites. These operators are given in Eqs. (23) and (24) 
of the companion paper [lj. In reference Q it was found that in the limit U/t — ► oo the deviations AD = AR C , _i, 
Ai? Cj+ i, AR St -x, and Ai? Si+ i generated by application of a general ./V-electron operator On onto any eigenstate of 
the spin a electron number operator © are restricted to the ranges of the inequalities (60)- (63) of the same reference. 
Obviously, given one of the four ranges of values defined by these four inequalities, the other three are dependent and 
follow from the relations given in Eq. (56) of Ref. @- 

Let us consider the deviations AR C _i, Ai? Cj+ i, Ai? Sj _i, and Ai? s+ i in the expectation values of the four operators 
given in Eqs. (23) and (24) of Ref. [l| in the case that both the corresponding initial state \4>i) and final state \(f>i') 
are energy eigenstates of the ID Hubbard model in the limit U/t—* oo. Note that it is assumed that the final energy 
eigenstate \<fri>) is contained in the state On\4>i) where On is the general A/"-electron operator under consideration. In 
the case of these initial and final states the above-mentioned deviations are such that, 

AR aJa = {4>v\R a ,i a \4>v) - {(pi\R a j a \(pi) \ 

= (i>i'(U/t)\ V\U/t)R aJa V(U/t)\MU/t)) 

- (Mu/t)\ vHu/t) R a>la V(U/t)\MU/t)) ; 

= (^{U/t)\M a ^J^(U/t)}-{^i{U/t)\M a , ff MU/t)}, (28) 

for all values of U/t. In order to derive this result we used Eq. (52) of the companion paper Q]. Also Eqs. (54)-(57) 
of the same paper were used in the evaluation of the equalities on the right-hand side of Eq. I|28|) . Equation (52) of 
Ref. |]| relates the energy eigenstates of the U/t — > oo Hubbard model to the corresponding energy eigenstates of the 
Hubbard model for any value of U/t. We note that the last quantity on the right-hand side of Eq. 12811 is nothing 
but the deviation 

AM a ^ a = (MV/t)\M a ,„Jil>t,(U/t)) - (MU/t)\M a ^ a \MU/t)) = M' attTa - M a , ffQ , (29) 

in the particular case when the initial state \tpi(U /t)) and final state \4>p {U/t)) are energy eigenstates of the finite-Z7 /t 
Hubbard model. The final energy eigenstate \tpi>(U/t)) is contained in the state Oj\j-\ipi(U/t)), where the rotated 
operator On is related to the operator On by Eq. I|26|) . The operator 

[On - N ] = [S, 6 N ] + \ [S, [S, On]} + ... , (30) 

can be expressed as a sum of A^'-electron operators such that AT' > N . Note that Eq. (|30|l contains the same 
information as Eq. (|26|l . 

The quantities M a , aa and M' a iJa of Eq. l(2"§|> are the eigenvalues of the energy eigenstates \ipi(U/t)) and (U / 1)) , 
respectively, relative to the holon (a — c) or spinon (a — s) number operator M Qi(Tq studied in the companion paper 
0. These two states are related to the above initial state \4>i) and final state | </>/') by the transformation given in 
Eq. (52) of the same paper. The equalities of Eq. I|28|) are valid for all possible 4^" initial energy eigenstates \4>i) of 
the ID Hubbard model in the limit U/t — > oo such that I = 1,2, ...,A Na and for final energy eigenstates \4>i>) which 
are contained in On\<Pi)- According to these equalities, the deviations in the eigenvalues of the four operators R a ,i a 
of Eqs. (23) and (24) of the companion paper .1] which count the number of electron doubly-occupied sites, empty 
sites, spin-down singly-occupied sites, and spin-up singly-occupied sites equal the deviations in the eigenvalues of the 
four operators M Q (Tq given in Eq. (22) of the same paper. These four operators count the number of —1/2 holons, 
+1/2 holons, —1/2 spinons, and +1/2 spinons, respectively. Note that according to Eqs. (52) and (54)-(57) of Ref. 
such an equality is valid because the initial and final states of the latter deviations are the energy eigenstates of 
the finite-?//* ID Hubbard model \ipi(U/t)) = V^(U/t)\4>i) and \ipi<(U/t)) = (U/t) | <j> v ), respectively. Here the final 
state \%l) V (U/t)) = V^(U/t)\4> v ) is contained in the state OaHVM^A)) = V\U /t) 6n\^i) ■ 

The electron double-occupation, no-occupation, spin-down single-occupation, and spin-up single-occupation eigen- 
value deviations (|28|l which result from application onto the state \4>i) of the above general operator On are restricted 
to the ranges of the inequalities (60)-(63) of Ref. @- Therefore, also the the — 1/2-holon-number, +l/2-holon-number, 
— 1/2-spinon- number, and +l/2-spinon-number eigenvalue deviations (|29|l which result from application onto the state 
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\ipi(U/t)} of the corresponding A/"-rotated-electron operator Oj^f are restricted to the ranges of the following inequali- 
ties, 



- J2 ^ AM c,-i/ 2 < J2 > ( 31 ) 

- ^ < AM c , +1/2 < ^ 7V+i,i a , (32) 
^ i-i,jV| e ,(.<AM Si _ 1/2 < ^ Si.,i.M e ,i., (33) 



; C ,Z S =±1 Z c ,Z„=±l 

and 



- I] i,i,4,i,<AM Si+1/2 < ^ h e ,- lm M e ,h- (34) 

According to Eq. (j29(l . the eigenvalues M a ,a- a and CT of these deviations refer to an initial state which is one of 
the 4^" energy eigenstates \ipi(U/t)) of the fmite-t//i Hubbard model and to a final state which is an energy eigenstate 
\4>v(U/t)) of the same model contained in the state Oj^\ipi(U/t)). However, since the 4^" states \ipi(U/t)) associated 
with the deviations of Eq. (|29() constitute a complete basis for the Hilbert space of the ID Hubbard model, the 
inequalities (|5T|) - (|54f> also apply to general deviations of the form, 

AM QiCTq = (^|M«,aJ^'> - &\M a , a Jij) , (35) 

provided that the arbitrary states \if>'} and \tp) are eigenstates of the spin a electron number operator ijSJ). 

It follows from the relations of Eq. (50) of the companion paper that out of the four inequalities l|3*T|) - (|S4Tl . only 
one is independent. By summation of the two inequalities H31|) and Il.'i2[l and of the two inequalities (|33|1 and l|34(l we 
find the following inequalities for the deviations in the total number of holons and spinons, 

-Af < AM C < Af; -Af < AM S < Af . (36) 

Note that the limiting values of the inequalities I136H are simply —Af and Af, where Af is the number of elementary 
rotated-electron operators of the expression of the operator Oj^. Finally, through the use of Eq. (28) of Ref. we 
find that the validity of these inequalities is equivalent to the inequalities 

-Af < AN C < Af ; -Af < AiVj < Af , (37) 

for the deviations AiV c = —AN^ in the number N c of c pseudoparticles and in the number of c pseudoparticle 
holes, respectively. 

We emphasize that the exact holon and spinon selection rules (|31|I - H34|I and l|3()l) refer to deviations generated by 
A/"-rotated-electron operators. We find below that in the case of deviations generated by the corresponding TV-electron 
operator the dominant holon and spinon microscopic physical processes lead to excited states that obey the inequalities 
(|3Hl - l|34ll and (|36|l . Such dominant microscopic physical processes control most of the A/"-electron spectral- weight. 
Interestingly, as in the U » t case, for U << t these rules are exact for the corresponding A/"-electron operator, the 
dominant processes becoming the only processes contributing to the A/"-electron spectral properties. On the other 
hand, for values of the on-site repulsion such that U « At these rules are not exact for A/"-electron operators but 
correspond to dominant processes that amount for more than 99% of the electronic spectral weight. For instance, 
in the case of one-electron creation operators, we find that the spectral weight generated by higher-order holon and 
spinon processes is extremely small. Such a weight corresponds to states with holon and spinon numbers outside the 
domains defined by inequalities l|31l) - (|34|l and (|36|l and is maximum for half filling and at U ps At where it amounts 
for about 0.25% of the total spectral weight. However, such a weight decreases for decreasing values of the electronic 
density. For instance, at quarter filling it is maximum also for U « At but amounts for 0.005% of the total spectral 
weight only. Moreover, for small electronic densities the dominant processes become the only processes for A/"-electron 
operators at any value of U ft. In the case of the one-electron creation operators the states generated by the dominant 
processes refer to a minimum amount of relative spectral weight for U ~ At. Nevertheless, such a minimum value 
corresponds for the dominant holon and spinon processes to 99.75%, 99.99%, and 100.00% of the total one-electron 
addition spectral weight for electronic densities n = 1/a, n = I /2a, and n « 1/a, respectively, as found in Sec. IV. 
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C. THE EFFECTIVE ELECTRONIC LATTICE 

The electron - rotated-electron unitary transformation performs a rotation in Hilbert space which maps electrons 
onto rotated-electrons such that rotated-electron double occupation is a good quantum number. As discussed below, 
the rotated-electron site distribution configurations that describe the energy eigenstates are independent on the value 
of the ratio U/t. In contrast, the electron site distribution configurations that describe these states are dependent on 
U. This is consistent with the U/t dependence of the electron double-occupation quantities studied in Ref. 

As discussed in the companion paper [lj , the electronic lattice remains invariant under the electron - rotated-electron 
unitary transformation. However, in order to distinguish the rotated-electron from the electronic site distribution 
configurations it is useful to introduce an effective electronic lattice. The rotated electrons occupy the sites of the 
effective electronic lattice, whereas the electrons occupy the sites of the real-space lattice. Such an effective electronic 
lattice has the same number of sites j = 1,2,3, ...,N a , lattice constant a, and length L = N a x a as the real-space 
lattice. 

In the companion paper [2j it is found that the rotated-electron site distribution configurations that describe the 
energy eigenstates include separated charge and spin sequences. The charge (and spin) sequences correspond to 
distribution configurations of the rotated-electron doubly occupied sites and empty sites (and spin-down and spin-up 
singly occupied sites). For the effective electronic lattice the number of —1/2 holons equals rotated-electron double 
occupation and thus plays the same role as electron double occupation in the real-space lattice. For finite values of 
U/t, the number of electron doubly-occupied sites, empty sites, spin-down singly-occupied sites, and spin-up singly- 
occupied sites of the real-space lattice electron site distribution configurations that describe the energy eigenstates, is 
not a good quantum number. In contrast, for the effective electronic lattice the number of rotated-electron doubly- 
occupied sites, empty sites, spin-down singly-occupied sites, and spin-up singly-occupied sites is a good quantum 
number for all values of U/t. This number equals the number of —1/2 holons, +1/2 holons, —1/2 spinons, and +1/2 
spinons, respectively, of any energy eigenstate. 

A property of major importance is that the effective electronic lattice site distribution configurations of the rotated 
electrons which describe the energy eigenstate \ipi(U/t)) of the relation (52)of the companion paper pj are the same 
as the corresponding real-space lattice site distribution configurations of the electrons which describe the state | </>;). 
Since this latter state is an energy eigenstate of the ID Hubbard model in the limit U/t — > oo, the rotated-electron 
site distribution configurations which describe the energy eigenstates are independent of the value of U/t and are 
precisely the same as the corresponding electron site distribution configurations of the U/t — > oo ID Hubbard model. 
In the companion paper 2j it is shown that this property is behind the fact that the c and a, v pseudoparticle band- 
momentum occupancy configurations which describe the energy eigenstates are also independent of the value of U/t, 
only the energy spectrum of such configurations are U /t dependent. The fact that the U/t — > oo selection rule ranges 
given in Eqs. (60)-(63) of Ref. Q are precisely the same as the finite-t//i ranges imposed by the inequalities i|31|l - ll34l) 
results from the U/t independence of the rotated-electron site distribution configurations which describe the energy 
eigenstates. 

Another important property is that for electronic densities smaller or equal to one, a ground state has zero rotated- 
electron double occupation. The same occurs for the ground state electron site distribution configurations of the 
real-space lattice in the limit U/t — > oo, the ground-state electron double occupation being in that limit given by 
D° = R® —1 = 0. This property is consistent with the above property that the rotated-electron site distribution 
configurations which describe the energy eigenstates are U/t independent. Note that creation of a rotated electron 
onto the ground state either transforms an empty site of the effective electronic lattice into a singly-occupied site or 
a singly-occupied site of the lattice into a doubly-occupied site. These two alternative transitions lead to —1/2 holon 
number deviations such that AM C _i/2 = and AM C _i « = 1, respectively. Therefore, the selection rule of Eq. H46|) 
leads to permitted final states such that the values of these deviations are AM Ci _!/ 2 = 0, 1. For the effective electronic 
lattice the operators M c ._i/2, M c + i/ 2 , M S| _i/2, and M Sj+1 / 2 count the number of rotated-electron doubly-occupied 
sites, empty sites, spin-down singly-occupied sites, and spin-up singly-occupied sites, respectively. Each site of the 
effective electronic lattice can either be doubly occupied by two rotated electrons of opposite spin projection, empty, 
singly occupied by a spin-down rotated electron, and singly occupied by an spin-up rotated electron. The result of 
Ref. 0] that the —1/2 holons are zero spin-singlet combinations of two electrons of opposite spin projection whereas 
the +1/2 holons are two-electronic hole quantum objects, is consistent with the —1/2 holons and the +1/2 holons 
corresponding to rotated-electron doubly-occupied and empty sites of the effective electronic lattice, respectively. 

In the companion paper 0] it is found that the effective electronic lattice introduced here is related to a set of 
effective pseudoparticle lattices. The concept of effective pseudoparticle lattice introduced in that reference arises 
from a local description of the pseudoparticles, alternative to the band-momentum q description extracted directly 
from the Bethe-ansatz solution in the companion paper pj. While in the limit U/t — > oo such a local description 
of the pseudoparticles refers to the real-space lattice, in the case of finite values of U/t the local pseudoparticle 
character refers to the effective electronic lattice, as discussed in Ref. 0. In that reference the local pseudoparticles 
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are described in terms of rotated-electron site distribution configurations. The concepts of a effective pseudoparticle 
elattice and of a local pseudoparticle introduced in the companion paper are used in the studies of the few-electron 
spectral properties of the model presented in Refs. [6fl l6lj. 

IV. THE DOMINANT HOLON AND SPINON MICROSCOPIC PHYSICAL PROCESSES FOR THE 

FEW-ELECTRON SPECTRAL PROPERTIES 

The commutation relations (60)- (62) of the companion paper imply that the six generators of the ?y-spin and 
spin SU(2) algebras given in Eqs. (2), (9), and (10) of the same paper and the momentum operator (J7J) have the same 
expressions both in terms of elementary electronic operators cj a and Cj, a and elementary rotated-electron operators 
Sj a and 5j tCr , respectively. Thus all these seven operators are both two-electron and two-rotated-electron operators. 
These two-electron operators are such that all operator terms of the second expression of Eq. Q26J1 vanish except the 
first term. Thus the operator (|30fl associated with such operators vanishes. This result suggests that application of 
an operator of the form (|5"U)) associated with a general ./V-electron operator onto a ground state leads to very little 
electronic spectral weight. The latter spectral weight vanishes in the particular case that the A/"-electron operator 
commutes with the electron - rotated-electron unitary operator, but is expected to be extremely small otherwise, at 
least for few-electron operators. The reason is that such general AA-electron and ./V-rotated-electron operators are 
products of the same elementary operators cj a , Cj jC and cj a , Cj >a , respectively, as the above two-electron operators. 
Application of any of the above seven two-electron operators and of the corresponding two-rotated-electron operators 
onto a ground state leads to the same final states. Other two-electron or few-electron operators have slightly different 
expressions in terms of the same elementary electronic operators. Thus it is expected that the application of any other 
few-electron operator onto a ground state leads to almost the same final states as application of the corresponding 
few-rotated-electron operator on the same ground state. This is equivalent to say that application of the operator l|3U|) 
associated with such operators onto a ground state leads to almost no spectral weight. This property is a direct result 
of the unitary and canonical character of the electron - rotated-electron transformation. This prediction is confirmed 
by the quantitative results obtained below. Indeed, we find that the dominant microscopic physical processes that 
control the A/"-electron spectral- weight distribution of the excitation Ojj\GS) are associated with the transformation 
laws of the A/"-rotated-electron operator of the second expression of Eq. <|26|) . Indeed, in the case Af — 1, over 99% 
percent of the spectral weight of Oj^\GS) is found to be contained in Ojj\GS). 

We start by studying the specific form that the exact holon and spinon selection rules associated with the inequalities 
(|31[I - <|34[1 take for transitions generated by application onto a ground state of a general A/"-rotated-electron operator. 
While these ground-state selection rules refer to A/"-rotated-electron operators, spectral functions of physical interest 
are of the form Hll(l - (|12(l and (|24|I - (|25JI and involve ./V-clcctron operators. Therefore, the main subject of this section is 
the study and discussion of the consequences of the A/"-rotated-electron ground-state selection rules on the A/"-electron 
spectral-weight properties. We find that the holon and spinon deviation number restrictions of the inequalities i|31|) - 
(f 341) correspond in the case of excitations generated by ./V-electron operators to the dominant holon and spinon 
microscopic physical processes that control the electronic spectral properties. 

In the quantitative numerical studies presented below we devote particular attention to the basic Af = 1 electron 
and rotated-electron operators. We provide evidence that the general results obtained for one-electron operators also 
apply to Af = 2 and other few-electron operators. 

A. THE GROUND-STATE CHARGE AND SPIN SELECTION RULES FOR ROTATED-ELECTRON 

OPERATORS 

The inequalities l|31|l - 134(1 provide the largest domains for permitted values of ±1/2 holon and ±1/2 spinon number 
deviations generated by application of ./V-rotated-electron operators O^f. According to Eq. (|lTj|) the number TV is 
given by Af = J2i i =±i ■M e , h ; where the numbers Af—i, — i, jV-i, +i, TV+i, -i, and A/+i,+i are defined below such 
an equation. The deviation permitted values defined by inequalities (|31[1 - <|34[1 are determined uniquely by the values 
of the four numbers Afi c j e specific to the operator Oj\f and do not depend on the values of U/t, electronic density 
n, and spin density m. Moreover, these inequalities provide the largest domains for permitted deviation values for 
all values of U/t, electronic densities < n < 1/a and 1/a < n < 2/ a and spin densities —n < m < n and 
— (2/a — n) < m < (2/a — n), respectively. However, the domains of such a permitted values can be smaller. Let 
M° ±1 y 2 be the numbers of ±1/2 holons (a — c) and ±1/2 spinons (a = s) of the initial state. For instance, in case 
of deviations which do not obey the following inequalities, 



12 



i,=±i i»=±i 

M°_ 1/a > £ S^-LMd.; Af° +1/a > £ 6i B ,i.Mio,i., (38) 
z. c ,z s =±i z c ,z 3 =±i 

the corresponding domains of deviation values are smaller and are contained in those defined by inequalities ({3*T|) - ({3~4*jl . 

For simplicity in this section we consider initial ground states corresponding to electronic densities < n < 1/a 
and spin densities < m < n. As discussed in the companion paper one can generate from states with electronic 
densities < n < 1/a and spin densities < m < n the corresponding towers of states with electronic densities and 
spin densities belonging to the above extended domains. This is reached by repetitive application of the off-diagonal 
generators of the 77-spin and spin algebras given in Eqs. (9) and (10) of Ref. Q, respectively, onto the former states. 

In the particular case when the initial state is a ground state \GS) with values of the density and spin density 
belonging to the above domains, the state Oj^\GS) belongs to an electron ensemble space with electron and rotated- 
electron numbers, 

N = iV T + N l = N° + AN ; 7V T = iV T ° + A7V T ; N t = Nf + AiV x . (39) 

Here iV°, N®, and N® are the ground-state initial electron-ensemble space electron and rotated-electron numbers and 
the deviations are given by the general expressions of Eq. 1)27(1 . (We recall that the spin a electron numbers remain 
invariant under the electron - rotated-electron unitary transformation and thus the number of spin a electrons equals 
the number of spin a rotated electrons.) 

Let us first disregard the selection rules associated with the inequality (|31|) . In this case the use of Eq. 1(27(1 
and Eqs. (49) and (50) of the companion paper reveals that the deviations leading to the set of final states 
whose c-pseudoparticle, —1/2 holon, and —1/2 spinon numbers are compatible with the state 0^/\GS) generated by 
application onto the ground state of a general ./V-rotated-electron operator are such that, 

AM c ,_ 1/2 =0,l,2,3,..,JVj+ £(-0^,-1=0,1,2,3,...,^. (40) 

ic=±l 

The corresponding possible values of the —1/2 spinon number deviations AM S) _im are determined by the value of 
the —1/2 holon number deviation AM Cj _!/ 2 and are given by, 

AAf s ,_ 1/2 = Hc)Mc,-i ~ AM c ,_ 1/2 = AN l - AM c ,_ 1/2 . (41) 
z c =±i 

For given initial values of the spin a electron and rotated-electron numbers, only one of the two deviations AM Cj _ X j 2 
and AM Si _jy 2 is independent. Thus Eq. (|4*0l reveals that the number of different CPHS ensemble spaces spanned by 
states whose electron and rotated-electron numbers are the same as the ones of the state Oj\r\GS) is iV° + AN± = N±, 
with Ni — > 00 in the present thermodynamic limit. This is confirmed by Eqs. 1(40(1 - (141(1 . 

According to the expressions given in Eq. (C24) of the companion paper pj the number of —1/2 holons vanishes 
in the case of a ground state and thus the inequality M" _ 1/ , 2 > J^i =±i-A/+i,i 8 given in Eq. 1(38(1 is not met except 
when __L 1 A/+i ) ; s = 0. This shortens the domain of the deviation values relative to the largest permitted domains 
given in inequalities ((31(l - l(34|) . When the initial state is the ground state the — 1/2-holon, +l/2-holon, — 1/2-spinon, 
and + 1/2-spinon number deviations are restricted to the ranges of the following inequalities, 

0<AM c ,_ 1/2 < £ (42) 
i,=±i 

£ (!cWi a ,i. < AM Ci+1/2 < £ (43) 
z c ,i s =±i z s =±i 



- £ ^_, s JV ic ,, s <AM Si _ 1/2 < ^H)JV ic ,-i, 

l a ,l, = ±l l c =±l 



(44) 
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and 

- Yl S lo,is^h,h < AM S|+1/2 < H)4,«' (45) 
; e ,;„=±i i c =±i 

Thus in the case of the —1/2 holon deviations AM Ci _x/2 only the following integer values are permitted, 

oc 

AAf Ci _i /2 = AZi Cf _i /2 +5^i/AJV c , I/ = 0,l,2,... ) M-i,i s . (46) 

f=i i a =±i 

Here we used Eq. (30) of the companion paper Q to express the deviation AAf c -1/2 m terms of the deviations in the 
numbers of —1/2 Yang holons and c : v pseudoparticles. We emphasize that the exact ground-state charge selection 
rule (|46|) also limits the value of the deviations ANj} in the number of c pseudoparticle holes. Indeed such a selection 
rule is equivalent to the following inequality, 

]T {lc)M B ,i. < AA^ 1 <M. (47) 

l c , i»=±l 

Note that this inequality is a particular case for an initial ground state of the general AN^ inequality l|37|l . 

In general the ground-state selection rule (|46|l also implies that the maximum permitted values of the deviations 
A£ C) _i/2 and J2™=i v AN C . v are max{ AL C: _ 1/2 } = Y,i s =±\N-\ t i s and u\ast{Y^ =1 v AN C ^ V } = J2i s =±iN-i,i 3 , 
respectively, provided that the condition max{(AL c _ 1 / 2 + Y^=i v AiV c , „)} = _ ±1 TV-i, i s is respected. However, 
we note that the value rn&x{X}^Li v AN CmV } — =±i-^— can be reached only if the inequalities (i) L° c +1 / 2 + 
S; 3 =±i-^+i,i 3 — Y^i a =±iN-i,i s and (ii) M° ±1 / 2 > 7V-i, T i +-V+i,±i are met. Inequality (i) follows from the 
requirement that generation of a c, v pseudoparticle involves combination of the v new created —1/2 holons with a 
number v of +1/2 holons 0. The number L° c +1 ^ 2 gives the number of these +1/2 holons pre-existing in the initial 
state and J2i =±i-^+i, h gives the number of +1/2 holons created by annihilation of rotated electrons singly occupying 
sites in the initial state. Inequality (ii) results from two requirements: First the maximum value of S w =i v AiV Ci v is 
reached when each new created spin-down and/or spin- up rotated electron combines with a spin- up and/or spin-down 
rotated electron pre-existing in rotated-electron singly occupied sites of the initial state. This requires that M s ° ±1 / 2 
should be larger than or equal to A/"_i, T i; Second at least A/+i, T i of spin +1/2 rotated electrons singly occupying 
sites of the initial state should be annihilated and give rise to empty sites associated with the new created +1/2 
holons. This requires the presence of the extra term TV+i, ±1 on the right-hand side of the above inequality (ii). 

For instance, for the half-filling 2S° = L® = ground state the values of the c, v pseudoparticle number deviations 
are restricted by the inequality Y^Li v AiV C! „ < min{^ ; __|_ 1 jV_i ) i s , J2i =±i-^+i,( s }- However, there might exist 
other selection rules which further restrict the values of these c, v pseudoparticle deviations. Recent numerical results 
[(57j reveal that in the case of the half-filling 2S® = lP c = initial ground state, the state generated by application of 
the J\f = 2 charge operator contains no c, 1 pseudoparticles. This suggests that there is another selection rule which 
states that generation of c, v pseudoparticles by new created rotated electrons results from combination of the —1/2 
holons generated by creation of these rotated electrons with pre-existing +1/2 Yang holons in the initial ground state. 
Since for the half-filling initial ground state we have that 2S® = LP C = 0, the numerical results of Ref. [67j suggest that 
in the case of such an initial ground state Y^Li v AiV Ci v — and thus Eq. (|46l) simplifies to AM C _jy 2 = AL C _ x y 2 = 
0, 1, 2, Y2i = ±i ■A/'-i, ( s ■ Moreover, if such a selection rule is valid for all ./V-rotated-electron operators, the maximum 
permitted value of the deviation Yl^Li v A-/V Cj v becomes max{^™ =1 v AN Ct „} = min{^ ; =±1 A/'-i, i s , LP, +1/2 }. This 
gives zero in the case of the half-filling initial ground state, in agreement with the numerical results of Ref. [67j • 

The number v equals the length of the ideal charge Takahashi's string excitations associated with the c, v pseu- 
doparticle branch [l| . We emphasize that an infinite number of the matrix elements between the initial ground state 
and excited states of the Lehmann representation of the A/"-rotated-electron spectral function of the operator Ojj van- 
ishes exactly as a consequence of the ground-state charge selection rule associated with Eq. I|46|). Moreover, we find 
later in this section that the exact A/"-rotated-electron rules correspond to the dominant holon and spinon microscopic 
physical processes generated by application of the corresponding TV-electron operator onto the ground state. For 
finite values oiU/t less than 1% of the TV-electron spectral weight does not correspond to the final states associated 
with deviation values obeying the selection rule l|46|) . Such a small amount of spectral weight corresponds to a few 
extra final states with deviation values outside the ranges given in Eq. J2BJI. However, it is found in Ref. [6l| that 
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the absolute value of the matrix elements of the spectral functions (|11|) - (|12|) and l|24l) - (|25l) decreases rapidly as the 
number of holon and spinon processes needed for generation of the corresponding excited states increases. It follows 
that also for the TV-electron operator an infinite number of possible final states give no measurable contribution to 
the A"-electron spectral weight. This property is related to the exact ground-state selection rule (fH))) for the states 
generated by the corresponding A/'-rotated-electron operator. 

In the case of A/'-rotated-electron operators, the absolute maximum value of the permitted values given in Eq. 
(Eht|i is given by Af. The value M C) _i/2 = Af is reached when the operator CV can be written as the product of Af 
creation rotated-electron operators. Analysis of Eq. (|4(j[) reveals that the general A/'-rotated-electron (and A/"-electron) 
excitation Oj^\GS) has quantum overlap (and Oj^\GS) has significant quantum overlap) only for energy eigenstates 
described by charge Takahashi's ideal string excitations of length v < J^i =±1^-1, ^ n the particular case when 
On can be expressed as a product of Af creation rotated-electron operators this gives v < Af. The number Af is 
thus the maximum absolute value for the length of such a charge ideal string excitation that can be generated by 
application onto the ground state of a general A/'-rotated-electron operator. In the case of the excitation generated 
by application onto the same state of the corresponding A/"-electron operator more than 99% of the spectral weight 
is exhausted by excited states containing string excitations of length v < _ ±1 Af-i, i a , as confirmed below for 
one-electron operators. This also gives v < Af if the A/"-electron operator can be written as a product of Af creation 
electronic operators. 

When the initial state is a ground state the permitted values of the —1/2 spinon deviation AM^.^ generated 
by application of the ground state of a A/'-rotated-electron operator are given by Eq. (|41|l . with AM C _i / 2 restricted 
to the values of Eq. This leads to a domain of values obeying the inequality 144|) . However, this exact spinon 

selection rule does not restrict the length v of the spin Takahashi's ideal string excitations which have quantum 
overlap with A/'-rotated-electron excitations Oj^\GS). This is because according to the expressions given in Eq. (C25) 
of the companion paper 0], the ground state occupancy of —1/2 spinons is finite and given by M s ° _ ± , 2 = iV° l = N®. 

Thus one can generate both —1/2 HL spinons and s, v pseudoparticles belonging to branches such that v > 1 by 
decomposition processes of s, 1 pseudoparticles which do not change the net number of —1/2 spinons. We note that 
the exact spinon selection rules just limit the value of the deviations of such a net number and do not limit the 
processes which conserve its value. 

The exact selection rule (|46[l is equivalent to the inequality (|47|l . This latter inequality reveals that the deviation 
in the number of c pseudoparticle holes generated by application onto a ground state of a rotated-electron operator 
is limited by an exact selection rule. Based on the expressions provided in the companion paper , one can express 
the number of c pseudoparticle holes and the number of s, 1 pseudoparticle holes in terms of the values of ry-spin S c 
and spin S s as follows, 

00 00 
JVj =2S c + 2j2vN c>1/ = L c + 2j2vN c ,„, (48) 

and 

oc oo 

N^ 1 = 2S. + 2^ f (u-l)N a , v = L a +2^ i (y-i)N a , v , (49) 
i/=i i/=i 

respectively. Based on the symmetries between the charge and spin sectors of the ID Hubbard model, one would 
expect that the values of the deviation in the number of s, 1 pseudoparticle holes generated by application onto a 
ground state of a rotated-electron operator should also be limited. Indeed, below we find numerical evidence that 
within the permitted final states which obey the selection rules (|46|) and (|47() , there is a sub class of states that describe 
over 94% of the few-rotated-electron spectral weight and whose deviation values in the number of s, 1 pseudoparticle 
holes are such that, 

- (lc-h)M e ,i B < AJVj'i <Af. (50) 

Note that the exact selection rule (|46|l limits the number of —1/2 Yang holons and c, v pseudoparticles generated 
by application onto a ground state of a AA-rotated-electron operator Oj^f. Also the states whose deviation values obey 
Eq. H50(l have the following restrictions in the number of —1/2 HL spinons and s, v pseudoparticles belonging to v > 1 
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branches 



oo 



AM s ,_ 1/2 = AL s ,_ 1/2 + ^(i/-l)AiV s , l/ = 0,l,2,..., ^ S lcJs M c j s . 



(51) 




Here the quantity M s< -1/2 is given by, 




M Sj _ 1/2 = L s , _i /2 + 5^(1/ - 1) iV a ,„ = M s , _ 1/a - ^ JV«, 



(52) 




Note that the value of this number vanishes for the ground state. Equations (|50|1 and (|51|) contain the same information. 

Let us consider the specific case of one rotated-electron removal and addition. In this case the states obeying the 
deviation value restrictions of Eqs. 1(50(1 and 1(51(1 are such that AN^ 1 = 1. Since the rotated-electron spectral weight 
distributions are independent of the value of U/t and as U/t — > 00 rotated electrons equal electrons, we can find 
the latter distributions by evaluating the electron spectral weight for U/t — > 00. By use of the method reported 
in Ref. [2^], we find that as the number of sites L increases the one s, 1 pseudoparticle hole excitations have for 
one-electron removal and addition and U/t — > 00 the relative weights provided in the Table below. We expect that 
for L — > 00 the value of the relative weights for one-electron removal and addition given in the Table for smaller 
systems are as U/t — ► 00 above 98% and 94%, respectively. The numbers provided in the Table seem to confirm 
that the deviation value restrictions of Eqs. (|50|l and 1(51(1 refer to processes that lead to a substantial part of the 
rotated-electron spectral weight associated with the permitted states whose deviation values obey the exact selection 
rule J5SJ. The U/t — > 00 results of Ref. [22j were derived by a scheme where the spin degrees of freedom of the 
ID Hubbard model were described by a Heisenberg isotropic spin model. Within such a scheme the one-electron 
excitations were in part simulated by a procedure involving a change in the number of sites of the Heisenberg chain. 
This procedure has similarities with the one used in Refs. [5(3, The deviation value restrictions of Eq. (|5 1|) do not 
correspond to an exact selection rule but are behind the result considered surprising in Ref. [2^ that for U/t — > 00, 
more than 97% and more than 99% of the spectral weight associated with creation of an electron and annihilation of 
an electron, respectively, is found on the spin branch of Faddeev and Takhtajan. This branch dispersion refers to the 
part of the one-rotated-electron spectral weight associated with the spin degrees of freedom. It is well known from 
the Bethe-ansatz solution that the holes in the c and s, 1 pseudoparticle bands play a key role in the energy spectrum 
of the charge and spin excitations of the ID Hubbard model, respectively. Note that the ground-state selection rules 
(146(1 and deviation value restrictions 1)51(1 can be expressed in terms of the deviation values of the numbers of c and 
s, 1 pseudoparticle holes through the inequalities 1(47(1 and 1(50(1 . respectively. When Ojsf is the one-rotated-electron 
creation or annihilation operator, according to the ground state exact selection rule 1(47(1 and the deviation value 
restrictions of Eq. 1(50(1 . the maximum number of holes created in the c pseudoparticle and s, 1 pseudoparticle bands 
respectively is one. Thus if the initial ground state has zero spin density the single s, 1 pseudoparticle hole indeed 
corresponds to the well known spin branch dispersion of Faddeev and Takhtajan 49]. 

Similar results are expected to hold for excitations generated by other few rotated-electron operators. The numbers 
provided in the table were obtained by numerical simulations and seem to confirm that, in contrast to Eq. ((46(1 . Eq. 
((51(1 does not correspond to an exact selection rule for deviations generated by rotated-electron operators. Note that 
the maximum value on the right-hand side of Eq. ((51(1 is the same as the corresponding AM 8> -1/2 maximum value 
of the general —1/2 spinon number deviation inequality 1(33(1 . 

B. DOMINANT PROCESSES FOR THE FEW-ELECTRON SPECTRAL PROPERTIES 

Here we confirm that the ground-state charge selection rule ((46(1 for deviations generated by ./V-rotated-electron 
operators also defines the dominant holon and spinon microscopic processes that control the spectral properties of the 
corresponding A/"-electron operator. The dominant processes correspond to 99.75% - 100.00% of the whole electronic 
spectral weight and are associated with states whose deviation values generated by application of a A/"-electron operator 
onto a ground state obey Eq. 1(46(1 . Within such dominant processes, there is a sub-group of simple holon and spinon 
processes that correspond to about 94.00% of the electronic spectral weight. These latter processes are associated 
with states whose deviation values generated by application of a ./V-electron operator onto a ground state obey both 
Eqs. id^t and (|5T)l . (We recall that the deviation value restrictions of Eq. (|5T)l are not exact selection rules for 
rotated-electron operators.) 
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N rotated-electron removal rotated-electron addition 



6 


0.998792 


0.977515 


8 


0.997486 


0.972141 


10 


0.996277 


0.968088 


12 


0.995178 


0.964847 


14 


0.994176 


0.962156 


16 


0.993258 


0.959862 


18 


0.992409 


0.957867 


20 


0.991622 


0.956105 


22 


0.990886 


0.954531 


24 


0.990196 


0.953109 



TABLE I: The relative weight of the one s, 1 pseudoparticle hole contributions in the U/t — > oo limit for finite size systems of 
N electrons. This also gives the same relative weight in the case of excitations generated by one-rotated-electron operators for 
any value of U/t. 

We start by considering the specific case of one-electron addition. We recall that according to the results of the 
companion paper I, the number ofholons M c = M Ci -i/ 2 +M Cj+1 / 2 and the number of spinons M s = M St _i/ 2 +M s ^ +1 / 2 
are such that M c = [N a — N c ] and M s = N c , respectively. For simplicity, let us assume that the initial ground state 
has zero spin density. In this case the spectral-weight distribution associated with creation of a spin-up electron has 
the same form as the one associated with creation of a spin-down electron. Here we consider the former case. The 
spin-up electron operator is related by Eq. 1(26(1 to the spin-up rotated electron. Thus in the present one-electron case 
the selection rule 1(4611 and the deviation value restrictions given in Eq. I|51|l refer to creation of a spin-up rotated 
electron. From the use of expressions (46), (47), (49), and (50) of the companion paper Q] for the electron and 
rotated-electron numbers in terms of the holon, spinon, and c pseudoparticle numbers, we find that in the case of 
creation of a spin-up rotated electron the following transitions are permitted by the selection rule 146JI and by the 
equivalent inequality (|47|> : 

(i) Transitions such that -AiV c = AiVj = -1, AM c _ 1/2 = 0, AM s ._ 1/2 = 0, AM C = -1, and AM S = 1. The 
minimal excitation energy for such transitions is zero. Within these general transitions, the transitions that also obey 
the restrictions f50"|l and (|5l"Jl are such that AN S , i = and AN£ 1 = 1. 

(ii) Transitions such that ~AN C = AN 1 ? = 1, AAf c ,_ 1/2 = 1, AM 8i _ 1/2 = -1, AM C = 1, and AM, = -1. The 
minimal excitation energy for such transitions is E u . Within these general transitions, the transitions that also obey 
the restrictions JSUJl and (|5T)l are such that AN S , i = — 1 and AN^ l = 1. 

In the case of creation of a spin-up rotated electron the selection rule l|46l) imposes that AM c _ 1 / 2 = 0, 1 and thus 
that AN' c l < 1. Note that the general transitions (i) and (ii) obey such restrictions. Within these general transitions, 
the transitions that also obey the restrictions 1(51(1 are such that [AM s _ 1 / 2 — AN S , i] = and thus that AN]} x < 1. 

The simplest non-permitted transitions involve creation of three holes in the c pseudoparticle band and of three 
holons: 

(iii) Transitions such that -AN C = AiVj = 3, AM C) _ 1/2 = 2, AM s> _ 1/2 = -2, AM C = 3, and AM S = -3. The 
minimal excitation energy for such transitions is 2E U . These transitions are not permitted by the rule 1(46(1 because 
AM C , _i/2 = 2 > 1. Within these general transitions, the transitions that obey the restrictions 1)50(1 and 151(1 are such 
that AN S , i = -2 and AN' S \ 1 = l. 

Within the permitted transitions of type (ii), the simplest transitions that do not obey the restrictions ((50(1 and 
((51(1 involve creation of three holes in the s, 1 pseudoparticle band: 

(ii') Transitions such that -AiV c = AN^ = 1, AM C) _ 1/2 = 1, AM 8i _ 1/2 = -1, AM C = 1, and AM S = -1 and 
with AN S , i = —2, AN]} 1 = 3. The minimal excitation energy for such transitions is E u . These transitions do not 
obey the restrictions (|5T|) because [AM Si _i/ 2 — AA^i] = AL s ^_i/ 2 = 1 > 0. 

In the case of half filling we often shift the ground-state zero-energy level to the middle of the Mott-Hubbard gap. 
In that case the above minimal excitation energies 0, E u , and 2E U become 0, Emh/2, and 3£ , M_f//2, respectively, 
where Emh = E u at n = 1. Below we call three-holon states the final states associated with the transitions (iii) 
because they involve creation of three holons. These states also involve creation of three holes in the c pseudoparticle 
band. Moreover, we call threes, 1-hole states the final states behind transitions (ii'), which involve creation of three 
holes in the s, 1 pseudoparticle band. In contrast, the final states associated with both the transitions (i) and (ii) are 
one-holon and one s, 1-hole states. 

Let us now find out what the relative weight of the states (iii) (or states (iii) and (ii')) is relative to the weight 
of the states (i), (ii), and (ii') (or states (i) and (ii)) when all these states are generated by application onto the 
ground state of the spin-up electron creation operator. To assess the importance of the three-holon final states (iii) 
and three-s, 1-hole states (ii'), first we turn to exact diagonalization of small chains. The small-chain results are 
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not expected to be a good approximation for the evaluation of the thermodynamic-limit weight distribution if we 
consider the weight associated with each specific final state. However, here we are mostly interested in the relative 
spectral-weight sum rules of the permitted states (i), (ii), and (ii') versus the forbidden states (iii). Moreover, we also 
want to know the relative weight of the states (ii'). Fortunately, we find below that in the case of such sum rules, 
which involve contributions from a whole class of states, the small-chain results provide values for the relative weights 
which agree up to 99% with the corresponding thermodynamic-limit values. 

The full electron addition and removal spectrum for six sites with six electrons (half filling) is shown in Fig^for a 
relatively large value of U —12 t. As mentioned above, in the case of half filling we define the Fermi level in the middle 
of the Mott-Hubbard gap. The Hubbard bands at ±Emh/2 and ±3Emh /2 are well separated, and the weights at 
w ±3Emh /2 are orders of magnitude smaller then the contribution of the main Hubbard bands, centered around 
±Emh/2- The states centered around 3Emh/2 are three-holon states of type (iii). As a result of the half-filling 
particle- hole symmetry, there is a corresponding structure for electron removal centered around —3Emh /2. In the 
large-f/t limit the latter structure is associated with creation of two real-space empty sites. For the general U/t case 
such a structure results from creation of two rotated-electron empty sites. We recall that the M summations of the 
few-electron spectral functions (|24|l - 1)25(1 refer to such a number of rotated-electron empty sites. 

In Fig. |21 we have plotted the contribution of different final states to the sum rule in the case of half filling. For that 
reason, we have followed adiabatically the weights (matrix elements) of different states as we reduced U/t (which in 
the present case we did for this relatively small system), and summed the weight over the particular family of states. 
As we can see, the contribution of the three-holon final states (iii) to the total sum rule is largest at intermediate 
values of U w 4t, and it does not exceed 0.25% in the total sum rule. For large values of U/t it decreases as (t/U) 4 , 
while for small values as (U/t) 4 . The three-s, 1-hole contribution from the states (iv) is also less than 0.25 %, and for 
small values of U/t it goes as (U/t) 2 . We find below that the relative spectral weight of the final states (iii) decreases 
with decreasing density. For instance, at quarter filling such a weight is 2% of that of half filling and vanishes in the 
limit of vanishing electronic density. Thus at quarter filling and U « At the final states (i), (ii), and (ii') correspond 
to 99.99% of the total spectral weight and the states (iii) to sw 0.005% of such a weight. On the other hand, we 
know from the values of the above table that the relative weight of the states (ii') increases with increasing values of 
U/t. As we confirm below, this is not so for the states (iii), that remain having very little spectral weight as L — > oo. 
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FIG. 1: Electron addition (uo > 0) and removal (u) < 0) spectrum for the half-filled six-site ring, with U — 12 1. Note the 
logarithmic scale for the weights. Not shown are the five-holon states, whose contributions are extremely small and energies 
are out of the shown energy window. 

While we couldn't do finite size analysis due to lack of larger system sizes for arbitrary value of U, for large U/t 
we have another method to calculate the contribution of the three-holon states (iii). When U/t is large, it is possible 
to achieve a more precise statement through the systematic t/U expansion of the electron - rotated-electron unitary 
operator [Tel ITTl [l2| . The electron - rotated-electron canonical transformation maps the Hamiltonian onto a rotated 
Hamiltonian. Such a transformation applies to all operators, as given in Eq. (|l(jfl . For instance, it also relates the 
elementary rotated electron operators to the elementary electron operators when calculating the spectral functions. 
So the creation operator according to the expansions ((15^1 and ljl"T|) can be decomposed into, 



U = 12 t 
n=1 



"].. a; -E u 



(53) 
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FIG. 2: The contribution of different states to the sum rule. Over 99% of the sum rule is exhausted by the one-holon and 
one-s, 1 pseudoparticle hole excitations corresponding to the above states (i) and (ii). For larger systems this remains true if 
we consider the above states (i), (ii), and (ii'), which are permitted states for the corresponding rotated-electron problem. 



where the index M in cl 



; ME, 



refers to the change of rotated-electron double occupation. Although in the present 



large U/t limit one has that E u w U, we use the general notation ME U for labelling each Hubbard band. The operator 

n hS +0{t/U) 



u j, ct; 



c'j a (l—nj t g-)+0(t/U) adds an electron to an unoccupied site. On the other hand, cj a . Eu 



3 ' ® 



adds an electron to a site already occupied by an electron of opposite spin, thus promoting the state to the subspace 
with one more doubly occupied site. In the U/t — ► oo limit these operators change the number of c pseudoparticles 
and holons accordingly to the deviation values of the above transitions (i) and (ii)-(ii'), respectively. 

The operator involved in the generation of the three-holon states (iii) is a . 2Eu ■ Its first non- vanishing term is 
proportional to (t/U) 2 and reads, 



L j,cr;2E u - jj2 



(54) 



where the rotated kinetic energy operators Tq and Te u are given in Eqs. (|2U[) and H21|) . respectively. After some 
algebra we find, 



t 2 
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V 5 f 




nj-i,(r(l-nj,a)nj+ 
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fl,cr4.ffS 




~ Hj—l, o)flj, cr^j'+l, cr 


-4 


cr4+l,ffS" 




-ftj-l,C7)^i,ff«j+l,C7 


-4 


<T C j-l,a C j 


+ 1, S n 3 


-1,ct"-j,ct(1-«j + 1.ct) 


+4 


-1, a^j+l, S 




-l,ff(l-«j,ff)«j+l,CT 


-4 




+i, s n i 


-1,ct"-j,ct(1-"j + 1.ct) 



(55) 

The operator only acts when all the j — 1, j, and J + 1 sites are singly occupied. The action of ct ^. 2B on sites 
j — 1, j, j + 1, apart from the t 2 /U 2 coefficient, gives 



TU) -> |dde) + |ded) 
III) -» -|dde) + |edd) 
lit) -> -|ded)-|edd) 
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while on any other configuration it will give zero. Here the indexes d and e stand for doubly-occupied sites and empty 
sites, respectively, and the | and | symbols refer to spin-up and spin-down singly-occupied sites, respectively. We can 
see that the operator ct ^. 2E removes three spins with total spin S=l/2 (for the S=3/2 combination | ][[) + | J,|| 

> + IUT>-o). 

Next, we calculate J A 21JHB (uj) du>, the total weight in the band centered around 3E u /2 w 3U/2 and of en- 
ergy width 12 t (each of the e and d sites contributes At to the bandwidth). It is given by the expectation value 



'j, a; -2B„ ( -j, cr; 2K 



| ^/>o)j where the is the spin-charge factorized wave function |17| . First we find that, 



/; 4 



(56) 



where the projector Pjj = hj t j + rij, j — 2nj f ensures the single occupation of site j. 
Replacing the factorized wave function, the partial sum rule reads, 



A 



:1< im (ij) dw = ( | - 2S S! - 2SiS 2 - 2S S 2 



(57) 



Let us denote by fj and fj the creation and annihilation operators, respectively, of a spin- less fermion at site j. 
The expectation value to find three neighboring spinless fermions is 



(nonin 2 ) B f 



(/oVo) (/oVr) (/ V 2 ) 
(/l/o) (AVr) (/i7 2 > 
(/l/o) (/l/i) (/l/ 2 > 
3 2nsin 2 (7rn) sin 2 (7m) sin(27rn) 



n sin 2 (27m) 
4^ 



(58) 



with the limiting behavior, 



(n nin 2 > sf 



47T b 



if n < 1 



135 

1 - 3(1 - n) if 1 - n < 1 



(59) 



Note that the spectral weigh of these states decreases rapidly away from half filling, i.e. at quarter filling (n — 1/2) 
it is about 2% of that at half filling. 

The expectation values in the thermodynamic limit of the Heisenberg model is |69l| , 



(SoS^Heis = ~-ln2« -0.443147 
(S S 2 ) H cis = i-41n2+^C(3)« 0.182039 



(60) 
(61) 



so that 



- - 2S Si - 2SiS 2 - 2S S 2 ) = 12 In 2 - -C(3) ~ 2.91. 

2 / Hcis ^ 



(62) 



The sum rule in the second upper Hubbard band defined in Sec. Ill is then, 



4 



.4 2UHB (w) dcu « 2.91 (hoh in2 } si 1— . (63) 



For the six-site finite-size cluster, the expectation value in Eq. (|62|l is (169 + 17\/l3)/78 s» 2.95, which is about 1% 
off from the thermodynamic-limit value. The asymptotic 2.95 t 4 /U 4 is shown in Fig. [3 as a dashed line. 
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Both our numerical results for general values of U/t and a small system and our analytical large U/t results for 
the thermodynamic limit confirm that the ground-state selection rule 146|) which limits the values of the deviations 
generated by rotated-electron operators defines the dominant holon and spinon microscopic processes that control the 
spectral- weight distribution of the corresponding electronic operators. While for U « t and U » t such a selection 
rule is exact also for these electronic operators, we find that the relative spectral weight of the permitted above states 
(i), (ii), and (ii') is minimum for U f=a At. This minimum value decreases with decreasing density. For half filling it 
is given by w 99.75%, whereas for quarter filling it reads ~ 99.99% and in the limit of vanishing density it becomes 
1=3 100.00%. The extremely small amount of missing spectral weight corresponds mainly to the forbidden three-holon 
states (iii). Higher order five-holon/five-c-hole states lead to nearly vanishing spectral weight. 

While the selection rule (|4(j|) is exact for rotated-electron operators, the states obeying the deviation value re- 
strictions of Eq. (|51|l correspond to rotated-electron processes that generate most spectral weight of the permitted 
transitions. Although the relative spectral weight of the three-s, 1-hole states (ii') increases with increasing the system 
length, we expect that for one-electron addition its maximum value is about 6% as L — > oo for half filling. However, 
higher order five-s, 1-hole states lead to nearly vanishing electronic spectral weight. Thus over 99% of the one-electron 
addition spectral weight corresponds to generation of one holon and one and three s, 1 pseudoparticle holes. 

A similar analysis can be performed for other few-electron operators. The deviation restrictions of Eq. (|46|) also 
amount to more than 99% of the spectral weight generated by application onto the ground state of A/"-electron 
operators such that Af > 1. Importantly, for all few-electron operators these deviation restrictions become exact 
both for U « t and U » t and refer to the dominant holon and spinon microscopic processes for U « 4t. The 
spectral weight associated with excited states whose deviation values are outside the ranges defined by Eq. i|46|) 
is for electronic operators in general maximum for half filling, yet it remains small for such an electronic density. 
As an example we consider the two-electron frequency dependent optical conductivity, which is directly related to 
the dynamical structure factor. A preliminary finite-energy study of such a Af = 2 electron problem was presented 
for the case of the ID Hubbard model in Ref. 63] . The ground-state selection rule of Eq. I|4t)l) implies that in 
that case the corresponding two- rotated-electron problem the permitted —1/2 holon number deviations are such 
that AM c _ 1 / 2 = 0, 1 and thus the permitted holon number deviations are AM C = 0, 2, respectively. The finite- 
energy absorption considered in that paper results from transitions associated with holon number deviations such 
that (AA/ C _!/2 = 0; AM C = 0) and (AM c _i/ 2 = 1; AM C = 2). Thus contributions from excited states with holon 
number deviations such that (AM c _i/j = 2; AM C = 4) are expected to correspond to less than 1% of the two-electron 
optical conductivity spectral weight. 

As for the one-electron problem, it is expected that such an extremely small amount of spectral weight is maximum 
yet very small at half filling. Such a prediction is confirmed by the results of Ref. |72j. There such a finite-energy 
absorption was investigated for the sine-Gordon model. The absorptions represented for coupling constant 1 = 0.9 
by a solid line and a dashed line in Fig. 1 of the same reference correspond to one-soliton/one-antisoliton and two- 
solition/two-antisoliton contributions, respectively. The latter absorption is very small, being multiplied by 100 in the 
figure. In the limit of (3 2 — > 1 the sine-Gordon model acquires a r^-spin SU(2) symmetry and describes the ID Hubbard 
model at half filling and small values of U/t. Moreover, in this limit the solitons and antisolitons become +1/2 holons 
and —1/2 holons, respectively. Therefore, the results of that reference are expected to be very similar to those of 
the half-filling Hubbard model. As /3 2 — ► 1 the above one-soliton/one-antisoliton and two-solition/two-antisoliton 
transitions correspond to (AA/ C +1 / 2 = 1; AM C = 2) and (AM Ci _i/2 = 2; AM C = 4) deviations, respectively. The 
transitions associated with the latter deviations are forbidden for the corresponding two-rotated-electron operator. 
Thus in the case of the two-electron operator the contributions from these transitions are expected to amount for 
less than 1% of the spectral weight. Such a prediction is confirmed by the smallness of the corresponding absorption, 
which is represented in Fig. 1 of Ref. jT^. Similar results hold for other two-electron operators. 

We emphasize that equivalent results also hold for one-electron removal. It is generally accepted that one-electron 
removal corresponds to final states with a single extra hole in the c pseudoparticle band. However, we find that this 
general expectation is only true for one-rotated-electron removal. In the case of one-electron removal less than 1% of 
the spectral weight goes again to final states with three extra c pseudoparticle holes. In the particular case of half 
filling such final states lead to the structure of Fig^which is centered around —3Emh/2. 

V. DISCUSSION AND CONCLUDING REMARKS 

In this section we summarize and discuss the results obtained in this paper. Moreover, we discuss the application 
of the concepts and non-perturbative many-electron tools introduced here and in the companion papers [l], to the 
accomplishment of a program for the study of few-electron spectral functions at finite values of excitation energy. 
Such an application program is fulfilled in Refs. |6Ct l6l|. 
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A. SUMMARY OF THE RESULTS 

The holon, spinon, and c pseudoparticle description introduced in the first paper Q of this series and further studied 
in this second and in the third papers, can be used in the evaluation of matrix elements between the ground state 
and excited states and in the related problem of the derivation of the line shape for finite-energy few-electron spectral 
functions |§3> EH1 • The results of the present paper about the dominant microscopic physical processes that amount 
for more than 99% of the spectral weight generated by application onto the ground state of few-electron operators 
is useful for the success of such a program. Our theory also describes the higher-order holon and spinon processes 
associated with the remaining less than 1% of electronic spectral weight. Furthermore, in this paper we used the 
relation of rotated electrons to the quantum objects whose occupancy configurations describe the energy eigenstates, 
in the generalization of the concepts of lower Hubbard band and upper Hubbard bands for all values of the on-site 
repulsion. While the lower Hubbard band refers to zero rotated-electron double occupation final states, the Mth 
upper Hubbard band corresponds to the spectral weight associated with the excited states of rotated-electron double 
occupation M. For large values of the ratio U/t this definition of lower Hubbard band and upper Hubbard bands 
coincides with the usual one. 

According to the inequalities given in Eq. , the maximum number of holons and spinons generated or annihilated 
by application of AT-rotated-electron operators onto any eigenstate of the spin a electron number operator is given 
by Af. In the particular case of excitations whose initial state is the ground state, the ranges provided in Eq. i|46|) 
refer to an exact selection rule for the deviation values generated by application onto the ground state of such Af- 
rotated-electron operators. The occurrence of this selection rule is a direct result of the relation of rotated electrons 
to the quantum objects whose occupancy configurations describe all energy eigenstates of the model. In this paper we 
found that the excited states associated with these permitted deviation values correspond to over 99% of the spectral 
weight generated by application onto the same state of the corresponding A/"-electron operators. This means that over 
99% of the spectral weight generated by application onto the ground state of a A/"-electron operator corresponds to 
application onto such a state of the first operator term of the second expression of Eq. (|26|) . Such an operator term 
is nothing but the A/"-rotated-electron operator associated with the A/'-electron operator under consideration. This 
interesting result reveals that application onto the ground sate of an operator of the form given in Eq. 1)31) [I leads to 
very little spectral weight. The value Af also determines the absolute maximum value v = Af of the quantum number 
v of the c, v pseudoparticle branches which can have finite occupancy in final states generated by A/"-rotated-electron 
operators. The number v is also the length of the ideal charge Takahashi's string excitations associated with the c, v 
pseudoparticle branch. Therefore, the A^-rotated-electron excitations only couple to energy eigenstates described by c 
Takahashi's ideal string charge excitations of length v such that v < Af. For the corresponding A/'-electron excitations, 
the excited states described by c Takahashi's ideal string charge excitations of length v such that v > Af amount to 
less than 1% of the electronic spectral weight. It follows that the dominant processes generated by application onto 
the ground state of a A/'-electron operator originate a number of 2^-holon composite c, v pseudoparticles and —1/2 
Yang holons and a number of c pseudoparticle holes whose values are within the ranges provided in Eqs. I|4()|) and 
(I47|l , respectively. We recall that for electronic densities n < 1/a there are no 2^-holon composite c, v pseudoparticles 
and no —1/2 Yang holons in the ground state. Higher order processes originating a number of these quantum objects 
outside the range defined by Eqs. ifTOjl and (|4"7|) correspond to less than 1% of the electronic spectral weight. These 
processes are generated by application onto the ground state of the operator given in Eq. l|3Ufl • 

Moreover, although the deviation value ranges provided in Eqs. (f 5011 and l|51|) are not associated with an exact 
selection rule for rotated-electron operators, within the dominant processes that include all values for these deviations, 
they refer to quite important processes for the electronic spectral weight. These equations refer to ranges for the 
numbers of generated or annihilated s, 1 pseudoparticle holes and for the numbers of generated 2^-spinon composite 
s, v pseudoparticles such that v > 1 and —1/2 HL spinons, respectively. We note that for spin densities m < n there 
are no 2^-spinon composite s,v pseudoparticles such that v > 1 and no —1/2 HL spinons in the ground state. For 
instance, in the case of one-electron addition considered in the previous section, Eq. (|50|l tells us that these important 
processes involve generation of one s, 1 pseudoparticle hole. We found that these processes lead to about 94% of the 
electronic spectral weight, whereas excited states involving the creation of three s, 1 pseudoparticle holes amount for 
more than 5% of such a spectral weight. We thus conclude that microscopic processes involving generation of more 
than three s, 1 pseudoparticle holes lead to almost no electronic spectral weight. Although in the case of the ideal spin 
Takahashi's string excitations associated with the s, v pseudoparticle branches the ranges of Eq. I|51[) are not exact, 
they correspond to a sub-class of excited states that amount to an important part of the electronic spectral weight. 

The occurrence of dominant holon and spinon microscopic processes for the electronic spectral weight and how fast 
the contribution of higher order processes vanishes is related to a property found in Ref. [61| that plays a central role 
in the evaluation of the line shape for finite-energy few-electron spectral functions. While the dominant holon and 
spinon microscopic processes associated with Eqs. (|46|l and (I47|l amount for more than 99% of the few-electron spectral 
weight, the remaining less than 1% of spectral weight is mostly associated with final states involving generation of 
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a few more c, v pseudoparticlcs, —1/2 Yang holons, and c pseudoparticle holes relative to the maximum value of 
the ranges provided in these equations. Also the contribution to the electronic spectral weight from excited states 
involving the generation of s, v pseudoparticles such that v > 1, —1/2 HL spinons, and s, 1 pseudoparticle holes 
decreases rapidly as the number of these quantum objects increases. It follows that excited states involving creation 
of a large number of the above-mentioned quantum objects lead to vanishing spectral weight. In particular, in the 
present thermodynamic limit the spectral weight of excited states whose generation involves an infinite number of 
quantum-object elementary processes vanishes exactly. Thus only processes involving generation of a finite number 
of quantum objects contribute to the few-electron spectral weight. This defines the Hilbert subspace of few-electron 
excitations. The pseudofermion description introduced in Ref. |60| refers to such a Hilbert subspace. The method 
for evaluation of matrix elements between the ground state and excited states introduced in Refs. [6fJ,|61| takes into 
account all possible final states generated by a final number of quantum-object processes, as discussed below. 



B. THE EVALUATION OF FINITE-ENERGY FEW-ELECTRON SPECTRAL FUNCTIONS 

The holon, spinon, and pseudoparticle description introduced in the companion paper and further studied here 
and in the companion paper is applied in Refs. jsjj, EH1 to the study of the few-electron spectral functions. The 
studies of these references provide the few-electron spectral function line shapes for all values of excitation energy. The 
above description is a necessary condition for the successful fulfilment of the program for evaluation of few-electron 
spectral functions at finite values of the excitation energy. Such a program is fulfilled in Refs -lfifiL l6lj| by application 
of the new concepts and paradigms introduced in this paper and in its companion papers [lj,|2j as follows: 

i) - The first step of the evaluation of a few-electron spectral function of the general form 1)11(1 - 1112(1 or 124(1 - 125(1 
is the expression of the corresponding few-electron operator as the expansion in terms of rotated-electron operators 
given by the second expression of Ref. 1(261) ; 

ii) - Next, one uses the results of the companion paper Q and express the rotated-electron operators in terms of 
local pseudoparticles operators, Yang-holon operators, and HL-spinon operators. Note that after use of the relation 
((26(1 . the problem of the evaluation of the matrix elements of the few-electron spectral functions ((12(1 and I25|) is 
equivalent to the computation of matrix elements involving the elementary creation and annihilation operators of 
these quantum objects ; 

iii) - However, rather than in terms of the interacting pseudoparticles, it is more appropriate to use in the evaluation 
of such matrix elements the related non-interacting operational pseudofermion description introduced in Ref. [60| . 
Fortunately, the few-electron spectral functions can be expressed as a convolution of the pseudofermion spectral 
functions corresponding to each pseudofermion branch with finite occupancy in the excited states \M,j) and \M,j) 
of the spectral- function expressions Hll(l - I(12|) and 1(24 (1 - 1(25 |) . respectively; 

iv) - Finally, the evaluation of the spectral function for each pseudofermion branch with finite occupancy for these 
excited states uses the method applied to U/t — > oo spin- less fermion operators in Ref. 22]. The energy spectra of 
the above excited states appearing in the spectral-function expressions 1(12(1 and ((25(1 is obtained by direct use of the 
Bethe-ansatz solution and 77-spin and spin symmetries. The state summations of these spectral- function expressions is 
simplified by the fact that the absolute value of the matrix elements vanishes rapidly as the number of pseudofermion 
processes involved in the generation from the ground state of the corresponding excited states increases. Such a rapid 
vanishing is an important property for the evaluation of the finite-energy few-electron spectral functions. It is directly 
related to the occurrence of dominant quantum-object microscopic processes. 

The use of the relation 1(26(1 in the evaluation of finite-energy few-electron spectral functions [(UEO! is a consequence 
of the breakthrough of the companion paper concerning the following two issues: First, the clarification of the 
relation of rotated electrons to the Bethe-ansatz quantum numbers that label the energy eigenstates; Second, the 
association of these quantum numbers with objects whose occupancy configurations describe the energy eig enstates. 
The occurrence of dominant processes found in this paper is also important for the studies of Refs. [60l l6l| . In 
addition to contributing to the understanding of the non-perturbative microscopic physical mechanisms that control 
the few-electron spectral properties, our results are related to the rapid vanishing of the matrix elements as the number 
of processes that generate the corresponding excited states increases, as discussed above. 
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C. FINAL DISCUSSION AND CONCLUDING REMARKS 



The complexity of the finite-energy one-electron and two-electron physics of the non-perturbative ID Hubbard 
model explains why except for U/t — ► oo investigations pol l2ll |22| , there have not been many previous studies about 
this interesting problem. In this paper we found that excited states generated from the ground state by microscopic 
processes leading to deviation values in the ranges defined by Eq. I|46l) lead to more than 99% of the few-electron 
spectral weight. Moreover, a careful comparison of the line shape of the one-electron spectral weight for the trivial 
U/t — > limit with the one obtained in Refs. [20H2H 12^ for the limit U/t — > 00, seems to indicate that most spectral 
weight is located in the vicinity of pseudoparticle branch lines. Such lines are generated by a set of final states where 
one pseudoparticle or pseudoparticle hole is created for all its available band-momentum values q and all remaining 
pseudoparticles or pseudoparticle holes are created at their Fermi points. In addition to these elementary processes, 
the line shape in the vicinity of these branch lines is for electronic and spin densities within the ranges < n < l/a 
and < m < n, respectively, generated by pseudoparticle - pseudoparticle hole processes in the c and s, 1 bands. 
The one-electron spectral-weight distribution obtained in Ref. [7lJ by numerical simulations refers to intermediate 
values of the ratio U/t such that U/t = 4 and displays charge and spin branch lines, in agreement with our prediction. 
Furthermore, this expectation is confirmed for both one-electron and few-electron spectral functions by the exact 
results of Ref. [6l|. 

Consideration of all possible branch lines constructed in this way for the dominant excited states whose deviation 
values obey the ranges defined by Eqs. I|46[l and (|51|) for a given few-electron operator is expected to describe the 
main spectral lines observed in real quasi-one-dimensional materials. This expectation is confirmed in the case of the 
organic conductor TTF-TCNQ by the preliminary application of our results presented in Ref. [E^. The separated 
one-electron charge and spin spectral lines observed in TTF-TCNQ by photoemission studies agree for the whole 
energy band width with the theoretical branch lines constructed by the above described procedure. These preliminary 
results confirm the interest for the study of the unusual finite-energy spectral properties of these materials of the 
concepts introduced in this paper. 

Elsewhere the powerful method introduced in Refs. |6d l6l| is used in a more detailed description of the line shape 
observed in TTF-TCNQ jfj^]. While the preliminary results of Ref. (5^ provide the form of the charge and spin 
branch lines, the detailed dependence on the excitation energy u> of the line shape in the vicinity of these lines is 
presented in Ref. [§2 ■ Fortunately, this more detailed study also agrees with the line shape observed in the real 
experiment. Elsewhere the method constructed in Refs. |6(J . |61| by means of the concepts and non-perturbative 
many-electron theoretical tools introduced here and in the companion papers 0, is applied to the study of the 
two-electron dynamical structure factor for all values of U ft [(jij. Moreover, a study of the phase diagram of a system 
of weakly coupled Hubbard chains which combines our ID results with a Renormalization Group scheme will be also 
implemented |65| . The results of these studies describe many of the anomalous spectral properties observed in real low- 
dimensional materials. In particular, they lead to spectral lines similar to the ones observed in finitc-energy/frequency 
experimental investigations of the one-electron spectral weight distribution or organic metals [29ll58| and two-electron 
dynamical structure factor of Mott-Hubbard insulators [2^| . Our theoretical predictions also describe and successfully 
explain the microscopic mechanisms behind the phase phase diagram observed in quasi-one-dimensional materials 
40]. 

Finally, the further application of the results obtained in this paper and in its companion paper [ij to the explicit 
evaluation of finite-energy few-electron spectral functions requires the introduction of the concepts of local pseudopar- 
ticle and effective pseudoparticle lattice. The introduction of these concepts is the main goal of the third and last 
paper of this series, Ref. 
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